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ABSTRACT 

Small angular scale (high £) studies of cosmic microwave background anisotropies 
require accurate knowledge of the statistical properties of extragalactic sources at cm- 
mm wavelengths. We have used a 30 GHz dual-beam receiver (OCRA-p) on the Toruh 
32-m telescope to measure the flux densities of 121 sources in VSA fields selected at 
15 GHz with the Ryle Telescope. We have detected 57 sources above a limiting flux 
density of 5 mJy, of which 31 sources have a flux density greater than 10 mJy, which is 
our effective completeness limit. From these measurements we derive a surface density 
of sources above 10 mJy at 30 GHz of 2.0 ± 0.4 per squ are degree. This is consistent 



with the surface density obtained by Mason et al. 
of sources selected at a much lower frequency 



( 2009 1 who observed a large sample 
1.4 GHz). We have also investigated 



the dependence of the spectral index distribution on flux density by comparing our 
results with those for sources above 1 Jy selected from the WMAP 22 GHz catalogue. 
We conclude that the proportion of steep spectrum sources increases with decre asing 
flux density, qualitatively consistent with the predictions of de Zotti et al. (2005). We 
find no evidence for an unexpected population of sources whose spectra rise towards 
high frequencies, which would affect our ability to interpret current high resolution 
CMB observations at 30 GHz and above. 
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1 INTRODUCTION 

Extragalactic radio sources are a serious contaminant in 
observations of primordial Cosmic Microwave Background 
(CMB) anisotropies on small angular scales (high I; e.g.|Tof-| 



|folatti et al.||2005 i. In order to estimate the CMB temper- 
ature power spectrum with high accuracy, detailed infor- 
mation about the foreground population is required. CMB 
experiments can be designed to work at any of a range of fre- 
quencies from ~20 to ~300 GHz, since the spectrum peaks 
at ~ 160 GHz. Extragalactic radio sources are problematic 
over most of this range, with different source populations, 
from active galaxies to star-forming galaxies, dominating the 
contamination in different wavebands. However, information 
on these contaminating populations is limited, since large- 
scale high-resolution surveys of the radio sky are restricted 
to lower frequencies, as for example in the NRAO VLA Sky 
Survey (NVSS; |Condon et "aI1p98| at 1.4 GHz and the 



Green Bank 4.85 GHz survey (GB 6; [Gregory et al.|p96 K 
Model-dependent predictions of the contaminating sources 
at the frequencies of CMB experiments may be made (e.g. 



Toffolatti et al. 1998 de Zotti et al. 2005), but these can 



be unreliable. Source populations can exhibit a wide range 
of spectral indices so that severe selection effects may arise 
when extrapolating to high frequencies from the limited set 
of surveys on which these predictions are based. Thus, for 
maximum reliability in their cosmological interpretation, it 
is essential that CMB observations are complemented by 
deep and carefully-selected high-frequency radio source sur- 
veys. 

The CMB power spectrum has been measured over a 
wide range of angular scales, from the early low-resolution 
(7° per beam) work of the Cosmic Background Explorer 
(COBE, |Smoot et al.|199"2| , through the intermediate scales 
probed by e.g. the Very Small Array (VSA, see e.g. Dick- 
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inson et al. 2004), Boomerang (see e.g. Jones et aL||2006 1 
and the Wilkinson Microwave Anisotropy Probe (WMAP, 
see e.g. Hinshaw et al.|200 3 2009), up to the high resolution 
measurements of e.g. the Cosmic Background Imager (CBI, 
see e.g. Sievers et al. 112009 1 and the Arcminute Cosmology 



Bolometer Receiver Array (ACBAR, sec e.g. Rcichardt et 
al.|2009 1 . The problem of contamination by radio sources in- 
creases at high multipole order, I, since the power spectrum 
of point sources increases as I 2 . Source subtraction strate- 
gies are therefore becoming ever more important. This is 
particularly so because of the results obtained with the CBI 
at 30 GHz, which indicate an excess of power at high £ com- 
pared to best fit cosmological models, the so-called 'CBI 
excess' ( Mason et al.|200 3 Rcadh ead et al.[ 2004; Sievers et 



al.|2009 1 . There has been much discussion of this result with 
some new observations not detecting a significant excess of 
power at i > 2000 (|Sharp et al.|[2009] |Friedman et al.|2009 



Reichardt et al.|2009 l, though most either do not cover the 
same range of multipoles or are at different observing fre- 
quencies. 

One possibility is that the excess arises from a pop- 
ulation of point sources, another is that that the power 
comes from the Sunyaev-Zel'dovich (SZ) effect associated 
with clusters of galaxies. Sievers et al. (20091 interpret the 



CBI excess as due to the SZ effect, provided that the value 
of ag, the normalization of the power spectrum of the den- 
sity fluctuations in the Universe on scales of 8/i _1 Mpc, is 
significantly higher than that calculated from WMAP data 



( Komatsu et al. 2009 ) . However, a value of as consistent with 



the WMAP results may still be allowed when the increased 
variance and skewness of the SZ effect statistics within small 



areas of the sky are taken into account ( Peel, Battye & Kay 
~ Shaw et aTp 009). The issue of point source contain- 



2009 



ination has been addressed by Sievers et al. ( 2009 1 using 



30 GHz measurements of NVSS sources in the CBI fields 



made with the GBT by Mason et al. (20091. Sievers et al 



( 2009 1 state that source contamination cannot account for 



the excess. The GBT 30 GHz sample was, however, selected 
at 1.4 GHz from the NVSS survey, so it is possible that 
a population of sources with spectra rising towards higher 
frequencies could have been missed. 

In order to look for the CBI excess, the VSA was up- 
graded and reconfigured in 2005 to provide both higher sen- 
sitivity and higher resolution, becoming the VSA 'Super- 
Extended Array'. With the issue of point source contami- 
nation in mind from the outset, the VSA system includes a 
separate pair of telescopes operating as an interferometer for 
source subtraction. As the Super-Extended fields lie within 
those observed during the project's earlier phases (see|Dick^| 



inson et al. 2004 and Grainge et al. 20031, information is 



therefore available on the stronger (~ 20 mjy) sources in 
the VSA target fields from the 33 GHz flux density mea- 
surements made using these telescopes ( Cleary et al.|2005 |. 
However, as discussed by Sievers et al. ( 2009 1, sources at the 



mjy level could combine to produce an excess of power at 
high multipoles. For this reason it is important to extend the 
work of |Cleary et al. (20051 to lower flux densities, however 



this could not easily be achieved using the VSA source sub- 
tractor. Ideally one would like to carry out a 30 GHz blind 
survey of sources down to flux densities of a few mjy but this 
is technically very challenging so an alternative approach is 
required. 



In this paper, we present the results from a collabo- 
ration between the VSA and the One Centimetre Receiver 
Array (OCRA, see e.g.|Browne et al.|2000||Lowe et al.|2007 



Lancaster et al.|2007 Pazderska et al.|2009| ) teams to study 



the population of radio sources that affect deep VSA mea- 
surements. The VSA group have surveyed the CMB target 



fields at 15 GHz using the Ryle Telescope (RT, see Wal- 
dram et al.|2003 i. We have then used the OCRA-p receiver, 



a 30 GHz dual-beam system, sensitive to a single linear po- 
larization, mounted on the 32-m Toruh telescope to measure 
the flux densities of the RT sources at 30 GHz, which is sim- 
ilar to the VSA's operating frequency of 33 GHz. 

Studies at 30 GHz of sources selected at 15 GHz also 
provide new information about the poorly-studied proper- 
ties of radio source populations at high radio frequencies. 
Understanding these populations at higher frequencies is vi- 
tal for upcoming work, not only for studies of CMB primor- 
dial anisotropics, but also planned cluster surveys utilising 
the Sunyaev-Zel'dovich effect. 



2 THE 15 GHZ SAMPLE FROM THE RT 
SURVEY 

2.1 The observations 

The five fields which have been observed by the VSA Super- 
extended Array, designated VSA 1L, 2G, 3L, 5E, 7E, are 
centred at RA( h m s ), Dec(° ' ") J2000: 00 16 09.1, +28 16 
40 (1L); 09 40 52.7, +31 46 21 (2G); 15 38 35.3, +41 40 
17 (3L); 03 02 57.0, +26 11 44 (5E); 12 32 21.8, +52 43 27 
(7E). These areas lie within those surveyed earlier by the RT 



for the work of Cleary et al. ( 2005 1, but the central regions, 



out to a radius of 1°, have now been scanned more deeply. 
This radius corresponds approximately to the first null in 
the relevant VSA primary beam. 

Since the FWHM of the primary beam of the RT is 
very small compared with that of the VSA (6 arcmin for 
the RT and 72 arcmin for the VSA), we have again used 
a rastering technique similar to that explained in |Waldram| 
( 2003 1 . The raster observations are combined to form 



et al. 



so-called raster maps which are used to identify peaks as 
possible sources. Each candidate is followed up with a short 
pointed observation, either to establish a reliable flux density 
or to eliminate it as a false detection. More details of how 
we have developed our method to reach greater sensitivity 



appear in Waldram et al. ( 2009 1, which describes the deeper 



9C 15 GHz survey. 

Each field has a sky area of 3.14 deg 2 , apart from 3L 
in which a very small fraction of the field (< 1 per cent), 
close to its low declination limit, was not covered by the sur- 
vey. The total area for all five fields amounts to 15.69 deg 2 or 
0.004778 sr. We have not assigned individual error estimates 
to the source postions and flux densities. The positions are 
derived from the raster maps and, by examining a sample 
of sources with independent measurements on overlapping 
maps, we conclude that the positional accuracy is better 
than 10 arcsec. The flux densities are from the pointed ob- 
servations and the uncertainty in these is dominated by the 
random uncertainty in the flux calibration, which we esti- 
mate to be ~ 5 per cent. 
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Figure 1. Plot of the ratio (S T /S P ) of raster flux density to 
pointed flux density versus S p . The estimated completeness limit 
is rs 7mjy. 
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Figure 2. The differential source count for the 15 GHz sample, 
showing the estimated completeness limit and the fitted 15 GHz 
count from the 9C survey. 



2.2 Completeness 

For the purpose of this paper an important question to con- 
sider is the completeness of the 15 GHz sample. The raster 
maps used to identify the candidate sources have a range of 
noise levels, from as low as < 1 mjy to as high as 1.3 or 
1.4 mjy at the edges of some of the fields. This means that, 
although we can be sure that every source in the catalogue is 
genuine, since each one has been followed up with a pointed 
observation, we are less sure of the number of undetected 
sources at any flux density level. We have approached this 
problem by a method described in detail in the deeper 9C 
paper (Waldra m et al.|2009[ ). Basically, we plot the ratio of 
the raster to the pointed flux density, (S T /S P ), versus S p (see 
Figure [l],. At the higher flux densities, where we can assume 
the sample to be complete, this ratio has a median value of 
0.86 (see the above paper for an explanation of why it is 
less than unity) but for the weaker sources the ratio rises 
systematically with decreasing flux density. This is because 
at these levels the peaks detected on the raster maps are 
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Figure 3. A cross-scan measurement of J0958+3224. First, a scan 
of the source is made in elevation, during which one beam crosses 
over the source. A Gaussian function is immediately fitted to this 
scan and an initial pointing correction is calculated and applied. 
A scan in azimuth is then performed, where both beams cross 
over the source. This is fitted using a double Gaussian function, 
and the final pointing correction is calculated. The initial pointing 
error is easily noticeable as there is a difference in heights of the 
beams in the two scans. The one second samples from these scans 
arc shown. 



preferentially those boosted by a positive local noise contri- 
bution. By investigating the progressive increase in the ratio 
with decreasing pointed flux density, we conclude that our 
completeness limit is approximately 7 mjy. This is consis- 
tent with the differential source count for the sample (see 
Figure [5| which shows a significant fall-off compared to the 
9C survey single power-law fit at the lower flux densities. 
There are 65 sources above the completeness limit and 56 
below. 

The completeness limit of the 15 GHz sample has im- 
plications for the 30 GHz list. On the basis that very few 
sources have a rising spectral index ^ 0.5 where S v oc 
v a , we conclude that we should have detected all sources 
within the VSA fields stronger than 10 mjy at 30 GHz. 
Hence our list of sources stronger than 10 mjy forms a com- 
plete statistical sample at 30 GHz. 



3 OCRA OBSERVATIONS AND 

CALIBRATION METHODOLOGY 

OCRA-p was used to observe the RT sample between Febru- 
ary and October 2007, with additional measurements made 
during 2008 and early 2009 to complete the sample and 
check for variability. The pointing positions RT, accurate 
to better than 10 arcsec, are more than adequate for the 
1.2 arcmin OCRA-p beam. Observations of discrete sources 
are made using two methods. Bright sources (S30 GHz <; 
lOOmJy) can be measured using cross-scan observations, 
where a scan across the source is made first in elevation and 
then in azimuth (see Figure [3| . This method is used both 
to measure the flux density of target sources (see e.g. Lowe 
et al. 2007) as well as to determine the required pointing 



corrections for the telescope (discussed further below). As 
the sources in the VSA fields are generally too weak for the 
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Figure 4. An example of an on-off observation on J1538+4105 
(27.6 ± 1.2 mjy). The first integration is taken when the source 
is tracked by the secondary beam (marked as (-) on the plot) for 
15 seconds (half of the standard integration time of 30 s). Next, 
three integrations are made, alternating between the primary and 
secondary beam, followed by a final integration with the source 
in the secondary beam of the same duration as the first. The last 
step is an integration at a background position 0.5 deg away in 
azimuth (marked as 'bg'). 



cross-scan technique, a second method was used whereby the 
respective two beams, the primary 'on' beam and the sec- 
ondary 'off' beam, are pointed at the source alternately with 
a switching cycle ~50 sec for a period ~5 min, thus mea- 
suring the source flux density relative to the background on 
each side of the source. This method is also used for observ- 



ing galaxy clusters with OCRA-p (see |Lancaster et al. 20071. 
The switching pattern for this on-off method has evolved 
during the observations and the final version is shown in 
Figure [4] 

The planetary nebula NGC 7027 was used as our pri- 
mary flux density calibrator, and was observed at least 
once per day during the observations of the VSA sources. 
NGC 7027 is the brightest planetary nebula in the radio sky. 
Hafez et all] p008t report a flux density of 5.39 ± 0.04 Jy 
at an epoch of 2003.0 at 33 GHz with a secular decrease of 
— 0.17±0.03 per cent per year. The flux density at the epoch 
of 2008.0 is therefore 5.34 ±0.04 Jy. We calculate a flux den- 
sity of 5.47 ± 0.04 Jy at 30 GHz, using a spectral index of 
-0.119 | |Hafez et al.|2008| |, and use this to calibrate all of the 
30 GHz flux measurements presented in this paper. [Zrjlstra, 
van Hoof fc Perley] (|2008 ) have also obtained accurate data 



for the flux density of NGC 7027 from which we obtain a 
value at 30 GHz of 5.37 ± .28 Jy at the epoc h of 2008.0, 
which is consistent with the Hafez et al.| ( |2008[ ) result. 

Pointing calibrations are carried out every 30-45 min- 
utes (or every 5-6 on-off measurements). The calibration 
sources used were J0015+3216 (for field 1L), J0958+3224 
(2G), J1521+4336 (3L), J0245+2405 (5E) and J1219+4829 
(7E). 



- the measured flux densities of the sources are found 
from two different fitting routines. For the azimuthal cross- 
scans of the calibration sources, two Gaussian peaks and a 
quadratic background were fitted, and the average of the 
two peak heights is used. 

The on-off measurements of the VSA sources are fitted 
by taking a symmetric double difference of the data. In this, 
the average of an 'on' measurement Son has the average of 
the second half of the preceding 'off' measurement S'off.prev 
and the first half of the following 'off' measurement Soff.next 
subtracted from it, i.e. 



Sdd = Jy ^2 0- 5 ("^n 



0.5(S' off , p 

rev + Soff.ncxt)) : 



(1) 



where the sum is over all N independent off-on-off sets in 
the data. As the difference between the measured fluxes is 
twice the source flux, an extra factor of 0.5 is introduced. 
The internal error on the measurement is calculated by 



— y 



°"s on 



+ 
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V 4 



■(2) 



Prior to this, any 1 second data points greater than 3 sigma 
away from the average - representing spikes in the data - 
are removed. 

A noise diode is used for secondary flux density cali- 
bration. The ratio of outputs obtained with cross-scans of 
NGC 7027 and the noise diode are averaged for a 24 hour pe- 
riod; this diode calibration is then applied to the individual 
measurements. 

The change of the telescope gain with elevation was 
corrected for using a fit to observations of NGC 7027 made 
over a wide range of elevations. As described by |Lowe et al.| 
(20071, atmospheric absorption is corrected for using mea- 



surements of the system temperature at the zenith and 60 
degrees from zenith made at regular intervals during the 
observations. A simple 'flat earth' model in which the at- 
mosphere is assumed to be a continuous slab with constant 
thickness and optical depth was used. 

A significant number of measurements had to be dis- 
carded due to adverse effects of poor weather. Each source 
was observed multiple times, with individual measurements 
with large fit errors removed, giving an average of 8.3 ob- 
servations per source. The final flux density for the source 
was calculated using the weighted mean and corre- 

sponding standard error on the mean. To ensure the maxi- 
mum reliability of the results, the observations were reduced 
separately at the three institutions involved in the OCRA 
collaboration, with the resulting values compared and dis- 
crepancies identified and resolved. The flux densities given 
in the next section are from one of these reductions; the 
other reductions give essentially the same values. 

The combined uncertainty from these factors is ~ 8 
per cent. The systematic error due to the error in the flux 
calibrator, NGC 7027, is negligible at 1 per cent. As such, 
the final error s on the 30 GHz flux densities are calculated 
by °" = V°ld~ 



+ (0.08S) 2 



4 DATA REDUCTION 

As the observations are carried out in two different ways - 
cross-scans across the calibrators and on-offs for the sources 



5 RESULTS 

The source names and J2000.0 positions at 15 GHz are listed 
in Table [I] along with their flux densities at 15 GHz from 
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Table 1. Sources in the VSA fields. The positions are from the RT 15 GHz observations. Sources marked with a 'v' in the note column 
are variable (see Section |6.1[ l, and those marked with an 'e' are extended in FIRST maps (see Section |6.4[ l. Notes on individual sources 
marked with a * are given in the Appendix. See Section [5] for details. 
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Table 1 — continued 
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dN 




the sample including selection effects; it does not represent 
the intrinsic source count distribution. 



S 30 (mJy) 

Figure 5. The binned numbers of the source flux densities 
at 30 GHz, including Poisson errors (•>/ max(dA r , 1)) (solid his- 
togram). The simple power law model is also shown (dot-dashed 
curve). Additionally, a simulated observational histogram (dot- 
ted histogram) shows how negative flux densities can occur due 
to Gaussian errors. 



the RT observations and at 30 GHz from the OCRA-p ob- 
servations. We find flux de nsities for the source s at 1.4 and 
4.8 GHz from the N VSS ( |Condon et ^ [1998 f] and GB6 
( Gregory et al. 1996| P| catalogues by searching for sources 
within 15 and 60 arcsec, respectively, of the 15 GHz source 
positions. Flux densities for all but a few sources were found 
in NVSS, although 50 sources were not present in GB6 prin- 
cipally due to the higher flux density limit of this survey. We 
calculate the two-point 1.4 to 30 GHz spectral index for the 
sources, defined such that S v tx v a . 



5.1 Distribution of flux densities 

We plot the distribution of 30 GHz flux densities in Figure 
[5] The negative flux densities for several of the sources listed 
in Table fT] are consistent with our error estimates. This can 
be illustrated using a simple model of the set of observed 
sources. We assume that the sources have been selected ran- 
domly from an intrinsic single power law model described 
by A^(5*3o) = 200 (Sso/lmJy) -1 over the whole survey area, 
with a sharp minimum flux density limit of 5*30 = 1.7 mJy, 
which is the median observational uncertainty. Figure [5] also 
shows the model and the average histogram from 100 simu- 
lated sets of values that would be observed if the model is 
correct and if the errors are Gaussian errors as listed in Ta- 
ble [l] The histo gram of our observations shows reasonable 
consistency with the model, i.e. the negative flux densities 
are consistent with random Gaussian errors. 

Note that the 1.7 mJy flux limit used here is well below 
the completeness limit of about 10 mJy (see Section [2.2[ ), so 
this power law model only provides a rough description of 



5.2 Source confusion 



for- 



It is possible that some of the 'off' measurements 
tuitously fell on previously uncatalogued 30 GHz sources 
rather than background sky. An upper limit to the level 
of confusion can be made as follows, by assuming that our 
sample is complete down to the median observational un- 
certainty of 1.7 mJy. Since the real completeness limit is 
about 10 mJy, this should provide a conservative upper 
limit. Based on the estimate of 31 GHz source counts from 
Mason et all |2009|, N(S S i) ~ 17 deg-^Sgi/lmJy)- - 8 , 



we expect about 0.003 sources/arcmin above 1.7 mJy. A 
compatible value is obtained from the much smaller 30 GHz 
survey by|Coble et al.| (|2007|). The OCRA-p beam width is 



1.2 arcmin, so the chance of a single 'off' measurement land- 
ing on an uncatalogued astrophysical object is about 0.003. 
For any single source, the flux density is estimated as an av- 
erage from 3-26 (mean: 8.3) individual measurements, each 
of which uses half the signal from each of two opposite 'off' 
positions (see section EJ. The individual measurements for a 
source are generally made at different position angles, with 
3.1 arcmin between 'on' and 'off' positions. We have 121 
sources with 1009 observations of 'off' positions. Hence, the 
expected number of 'off' flux densities which are contami- 
nated by uncatalogued objects is about 3.5. The chance that 
more than one of four contaminating 'off' flux densities occur 
for the same source is « 4.8 per cent [^] assuming that clus- 
tering effects are negligible. Hence, it would be statistically 
expected that about four sources each have one contami- 
nated 'off' flux density measurement, and no sources have 
multiple contaminations. 

For any of the sources affected by one contaminated 
'off' flux density measurement, the averaging procedure for 
the multiple measurements is likely to make the effect of 
the contamination small, since on average 8.3 measurements 
are made, and a single 'off' flux only contributes half of its 
value to a given measurement. Hence, on statistical grounds 
we expect that confusion is negligible for our sample as a 
whole, though it is likely to have a small effect on a few 
percent of our sources. 

In practice we have direct knowledge of the confusing 
fields from the RT 15 GHz observations themselves and from 



the N VSS ( |Condon et al.||1998[ > and FIRST ( |White et al. 
19971 catalogues at 1.4 GHz. In a few instances, sources in 



the RT 15 GHz list are either two parts of the same source 
or unrelated sources sufficiently close to each other in the 
sky as to be confused in the present OCRA observations; 
these are identified by brackets in Table [I] We have exam- 
ined the NVSS images of all sources and, where available, 
FIRST images. In the following cases, we conclude that our 
measured 30 GHz flux densities should be treated with cir- 
cumspection: J0014+2852a/b, J1535+4142/1535+4143a/b, 
J0301+2541/42 and J1238+5249a/b. More details on these 
particular sources are given in the Appendix. 



Accessed via http : //www. cv.nrao . edu/nvss/ 
2 Accessed via http://heasarc.gsfc.nasa.gov/W3Browse/all/ 3 The calculation is equivalent to that of the Birthday Problem: 

|gb6.html| p = 1 - [(1 - 1/121) x (1 - 2/121) x (1 - 3/121)] »s 0.049. 
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Table 2. Comparison of the 15 and 30 GHz flux densities from the RT and OCRA-p (in mjy) from common sources measured with 
the VLA at 1.4, 4.8, 15.2 , 22 and 43 GHz by |Bolton et al.| | |2004) ; |Waldram et al. 1(2007) and those measured by the source subtractor 
jCleary et al.|2005| [20l38] ). 









VLA 






RT 


OCRA 


VSA 




Name 


Si .4 


Sa.s 


Sl5.2 


S22 


S43 




S30 


S33 


Notes 


J0013+2834 


32.7 


33.1 


34.6 


36.6 


30.5 


23.8 ± 1.2 


14.5 ± 2.0 


36.8 ± 1.8 


Variable 


J0014+2815 


80.4 


60.1 


45.5 


37.7 


23.8 


45.1 ± 2.3 


27.6 ± 2.8 


34.8 ± 2.6 


Variable (FigJHJl 


J0019+2817 


25.9 


23.8 


17.4 


27.0 


34.9 


32.6 ± 1.6 


26.0 ± 3.5 


23.9 ± 1.7 


J0303+2531 












55.1 ± 2.8 


45.6 ± 5.2 


30.0 ± 2.1 


J0937+3206 


108.9 


53.9 


58.4 


58.8 


41.0 


49.3 ± 2.5 


61.8 ± 7.9 






J0944+3115 












24.3 ± 1.2 


32.1 ± 3.1 


25.3 ± 2.6 




J1235+5228 












47.0 ± 2.4 


27.3 ± 2.6 


27.9 ± 1.8 




J1539+4123 












31.8 ± 1.6 


36.2 ± 3.3 


33.7 ± 2.4 




J1539+4217 


53.3 


40.0 


34.2 


37.0 


26.3 


31.7 ± 1.6 


32.1 ± 3.0 


34.0 ± 1.6 




J1538+4225 


42.0 


40.5 


41.6 


42.2 


29.0 


31.9 ± 1.6 


27.5 ± 3.4 


41.0 ± 1.5 


Variable 


J1540+4138 


16.0 


30.9 


34.2 


23.0 


9.1 


45.6 ± 2.3 


25.3 ± 3.3 






J1541+4114 


65.0 


38.0 


30.5 


27.5 


19.2 


37.8 ± 1.9 


15.9 ± 2.3 


32.5 ± 2.3 


Variable 




400 500 600 
JD - 2454000 



900 



Figure 6. The flux density measured for J0019+2817 as a func- 
tion of time. The source is clearly variable. 
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Figure 7. Comparison of the OCRA-p flux densities against 
those from the VSA source subtractor. 



6 DISCUSSION 

6.1 Comparison with previous measurements and 
source variability 

The stronger sources (> lOmJy at 15 GHz) within the 33 
extended VSA fields were observed during the initial VSA 
campaign at 33 GHz using the VSA source subtractor, with 
those greater than 20 mjy at 33 GHz used for direct sub- 
traction from the CMB data ( Dickinson et al.|2004 l. A sub- 
sample of these was published by Cleary et al. (2005 20081. 



Table [2] lists the flux densities measured by the source sub- 
tractor for the sources in Cleary et al. ( 2005 2008 1 that also 



lie within the super-extended array fields, as well as mea- 
surements of those sources made at a range of frequencies 
by |Bolton et aL] (|2004|> and |Waldram et al.| (|2007[) using 



the Very Large Array, and compares them to the measure- 
ments by the RT and OCRA. A number of the sources show 
evidence of variability. 

We note that OCRA-p is only sensitive to a single lin- 
ear polarization, which could result in discrepancies between 
measurements from the VSA, VLA and OCRA for polar- 
ized sources; however these are likely to be small since the 
WMAP sources selected at 22 GHz (see Jackson et al. in 



preparation) exhibit average linear polarizations of ~3 per 
cent. The sources observed with the higher resolution in- 
terferometers (RT and VLA) may also be slightly resolved. 
Nevertheless we identify at least four sources (J0013+2834; 
J1538+4225; J1541+4114; J0019+2817) as likely to be in- 
trinsically variable. One of these sources - J0019+2817 - 
shows clear variability from OCRA data alone (see Figure 
[6|, although the mean flux density from the measurements 
agrees well with the source subtractor measurement. 

A comparison of flux densities for 42 sources observed 
both with the VSA source subtractor and OCRA-p is shown 
in Figure [7] The flux density scales are consistent, however 
10 out of the 42 sources show variations greater than 2er 
(combined error), some of which show variations greater 
than a factor of 2. These are denoted with a V in Table [l] 



6.2 Source spectra 

Figure [8] shows the '2-colour' diagram for sources detected 
at all four frequencies (1.4, 4.8, 15, 30 GHz). We find few 
sources with peaked spectra (top left quadrant of Figure [sj , 
but a larger number of sources with spectra rising towards 
high frequencies (right hand quadrants of Figure [SJ> . 
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Figure 8. a*'| vs. a^® for all sources in this sample with a 30 GHz 
flux density greater than 10 mjy, and known flux densities at 
all four frequencies. The diagrams in the corners schematically 
illustrate the spectral behaviour of the sources in each quadrant. 
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Figure 9. The distribution of spectral indices between 1.4 and 
30 GHz (dj° 4 ) for 60 sources in this sample with a 15 GHz flux 
density greater than 7 mjy. 3cr detections at 30 GHz are shown 
in the light grey histogram; the dark grey histogram contains the 
3<r upper limits on the spectral indices. 



Figure [9] displays the distribution of spectral indices 
S oc u a between 1.4 and 30 GHz for 60 sources above the 
completeness limit of 7 mjy at 15 GHz (5 out of a total 
of 65 sources are excluded due to their complex morphol- 
ogy - 0258+2530a/b, 0938+3140 and 1236+5305/6). There 
is clear evidence of bimodality, with a split around a spec- 
tral index of -0.5, dividing the sources into groups of sources 
with 'flat' and 'steep' spectra. 

The distribution is markedly different from the 1.4 to 
30 GHz spectral index distribution presented by |Mason et| 



al. (20091. However, their source sample was selected at 



1.4 GHz rather than 15 GHz as in this paper. The distribu- 
tion in Mason et al. ( 2009 1 is not bimodal, but has a single 



peak around —1.0, implying that their sample is dominated 
by the steep spectrum sources. Our sample contains a larger 
fraction of flat and rising spectrum sources. This is not sur- 
prising given the difference of a factor ten in the sample 
selection frequency. 

We return to the issue of spectral indices in different 
samples in sections 6.3 and 6.4. 
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Figure 10. The spectral index distribution of 22 GHz-selected 
WMAP sources between 1.4 and 33 GHz. The majority of sources 
have flat spectra (a > —0.5). 



6.3 The 1.4-30 GHz spectral index distribution vs 
flux density 

We are seeking to improve our knowledge of the source popu- 
lation at 30 GHz down to mjy flux densities and beyond (by 
extrapolation). In order to extrapolate one needs to know if 
there is a dependence of the spectral properties of the popu- 
lation on flux density at 30 GHz. We have therefore used the 



WMAP 5-year point source catalogue (Wright et al. 20091 



to select a complete strong source sample consisting of all 
sources with a 22 GHz flux above 1 Jy. We have then cross- 
identified these with the NVSS catalogue and calculated the 
spectral index distribution between 1.4 and 33 GHz. We 
chose 22 GHz as the selection frequency to be as close as pos- 
sible to the 15 GHz selection frequency used in the present 
paper. 

In this analysis, we exclude sources where there are 
no obvious NVSS matches, sources known to be extended 
on scales greater than ~ 30 arcmin (WMAP J0322- 



3711/Fornax A; WMAP J1633+8226/NGC 6251) or con- 
fused in WMAP (WMAP J0223+4303/3C66 A and B). 
Sources are matched to single NVSS sources, except for 
WMAP J0108+1319/3C33 where four close NVSS sources 
are combined. In total, 211 sources are used in the analysis. 
Figure [10] shows the spectral index distribution, which has 
a single peak centered on zero. All but 16 sources are classi- 
fied as flat spectrum, i.e. have a > —0.5. As expected, this 
is also very different to the distribution found by |Mason et| 



al. (20091, which is dominated by steep spectrum sources. 



There is no evidence for a significant population of steep 
spectrum sources in the WMAP spectral index distribution. 
This contrasts with the VSA source distribution shown in 
Figure[9]where there is evidence for both flat and steep spec- 
trum populations. This difference could be ascribed to the 
selection frequencies (22 GHz cf. 15 GHz), however their 
closeness suggests that this should not be a major issue. 
Another consideration is that although we are selecting us- 
ing the 22 GHz flux densities, in reality the WMAP point 
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Figure 11. The distribution of spectral indices between 1.4 and 
30 GHz (ofJ° 4 ) for 29 sources in this sample with a 30 GHz flux 
density greater than 10 mjy. The darker boxes represent sources 
with flux density between 10 and 20 mjy at 30 GHz; the lighter 
boxes represent the sources with greater than 20 mjy at 30 GHz. 
The distribution of the stronger sources closely resembles that of 
the WMAP sample. 



source catalogue is obtained using all of the bands (Wright 
ct al. 2009). This should also not be an issue as the sources 
with 22 GHz flux density greater than 1 Jy are detected 
by WMAP at high significance. A final possibility, which 
we think is most likely, is that the dissimilar distributions 
are due to the difference in the flux densities of the sources 
in the samples. The population of steep spectrum sources 
increases as the flux densities decrease. 



Such an effect is also seen in Waldram et al. (20091, 



where they find a significant change in spectral index dis- 
tribution with flux density amongst their 15 GHz selected 
sources. We note, however, that there is some overlap in 
sources between their samples and the one studied here. 

The bimodality of the spectral index distribution re- 
mains the same when the sources in this sample are se- 
lected at 30 GHz. Figure [Tl] shows the spectral index dis- 
tribution of the VSA sources selected to be complete at 
30 GHz (S30 > 10 mjy), and Figure [l2"| shows the spec- 
tral indices of WMAP sources selected to be greater than 
1 Jy at 33 GHz. Splitting the VSA source sample complete 
at 30 GHz into high (>20 mjy) and low (<20 mjy) samples 
shows that the higher flux density sources have a spectral 
index distribution closer to the WMAP distribution; the 
lower flux density sources become increasingly more steep 
spectrum. A Kolmogorov-Smirnov comparison of the VSA 
(S30 > 10 mjy) and WMAP distributions shows that the 
samples are different at the 95 per cent confidence level. The 



semi-empirical model described by de Zotti et al. ( 2005 I pre- 



dicts a cross-over between the dominance of flat and steep 
spectrum sources at about 30 mjy (their Figure 14). Our 
observations display a sharper cross-over than the models 
appear to suggest, and at the slightly lower flux density of 
20 mjy, but are otherwise consistent with the predictions by 
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Figure 12. The spectral index distribution between 1.4 and 
33 GHz of WMAP sources selected to be stronger than 1 Jan- 
sky at 33 GHz. The distribution is essentially the same as for 
the 22 GHz selected sample, but with slightly fewer sources with 
spectral indices between and -1. 



6.4 Morphological properties 

Three of the VSA fields are in the region covered by the 
FIRST survey ( White et al.| |l997), and thus 5 arcsec reso- 
lution 1.4 GHz maps are available. From a visual inspection 
of these sources, with contours plotted down to 4 per cent of 
the peak, we find that 41 of the 68 independent sources (60 
per cent) show extended structure; these are marked with 
an 'e' in Table [l] Of these sources, the majority have steep 
spectra (aT° 4 < -0.5). 

Sixty two sources in the sample of WMAP sources de- 
scribed in the previous section also have FIRST maps avail- 
able. Again we made a visual inspection of maps with the 
contours drawn to 4 per cent of the peak flux density. We find 
that only 21 per cent of these sources are extended, a factor 
three fewer than in the VSA sample. This difference in the 
proportion of extended sources between samples containing 
weak and strong sources reinforces our conclusion, based on 
spectral indices alone, that the mixture of source types is 
different in the two samples. It is important to consider if 
this could arise from differences in selection other than the 
flux density limit. There are several factors which might bias 
the content of the samples, often in opposite directions: 



de Zotti et al. (2005). 



• The VSA sources are selected at a lower frequency and 
therefore there will be a small bias towards more steep spec- 
trum extended sources compared to WMAP. 

• The resolutions of the surveys are very different; the RT 
has a beam of ~30 arcsec at 15 GHz, whilst WMAP has a 
beam of 0.88° at 22 GHz. The RT observations will discrim- 
inate against sources with structures J; 30 arcsec i.e. against 
very extended steep spectrum sources because emission will 
be resolved out and the sources will drop below the selection 
limit. 

• Some VSA sources are weak and thus in several of the 
FIRST images extended structure would be more difficult 
to detect than in the much stronger WMAP sources. 

• We do not know the redshift distribution of the VSA 
sources but being weaker they are likely to be at higher 
redshift on average than the WMAP sources. This could 
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favour detection of extended structure in the WMAP sources 
since lower redshift sources will have larger angular sizes and 
the K-correction will favour higher surface brightness. 

On balance we think the known selection effects are 
slightly biased towards finding more extended WMAP 
sources. We conclude that the population mix really does 
depend on flux density. 

Is the result unexpected? The magnitude of the effect 
may be larger than expected but not the sign. There is evi- 
dence in the data presented in Orr & Browne ( 1982[ ) that at 
5 GHz the proportion of flat spectrum (and therefore com- 
pact) sources decreases with flux density and they show that 
this is a robust prediction of relativistic beaming models. A 
clearer observational demonstration of the change of pop- 



completeness limit, we derive a surface density of sources 
at 30 GHz of 2.0 ± 0.4 per square degree. This is consis- 
tent with the value obtained by Mason et al. (2009), who 



ulation with flux density can be found in Condon ( 1984 1 



Condon's Figures 8 and 9 show that the compact flat spec- 
trum source counts at 5 GHz fall off faster with decreasing 
flux density than the counts of extended sources. This trend 
probably does not continue down to //Jy levels since results 



on the Hubble Deep Field (e.g. Muxlow et al. 2005) indi- 



cate that compact AGN still form about 25 per cent of the 
population which at these flux density levels is dominated 
by star-forming galaxies. For studies of the CMB, the exact 
population mix is not particularly relevant for total intensity 
observations, however it may well be important for polarized 
observations since extended AGNs are likely to have higher 
linear polarizations than their compact counterparts. 

In summary, the above results indicate a marked de- 
crease in the proportion of flat-spectrum sources as a func- 
tion of decreasing flux density. Such a change in population 
mix may be important for the accurate prediction of polar- 
ized foregrounds for CMB experiments. 



6.5 Comparison of estimates of source surface 
densities at 30 GHz 

There are 31 sources in the five VSA fields detected above 
10 mjy at 30 GHz. These effectively comprise a complete 
sample as the sources were selected from a deeper survey 
at the nearby frequency of 15 GHz. The surface density of 
sources is thus 2 ± 0.4 sources per square degree (using the 
Poissonian error), which can be compared with the extrapo- 
lated value of 2.65 ± 0. 1 in |Mason et al] §2009\ . |Coble et al.] 
( 2007 1 also give an estimate of the surface density of 30 GHz 
sources in the field, which gives 2.2^i 2 sources per square 
degree; they also point out that the density of sources in 
clusters of galaxies is significantly higher. 



7 CONCLUSION 

In order to aid the subtraction of individual sources from the 
VSA fields observed at ~30 GHz, and to obtain a statistical 
estimate of the surface density of sources at 30 GHz, we 
have observed a sample of 121 sources using the OCRA- 
p receiver on the Toruh 32 m telescope; the sample was 
selected at 15 GHz with the RT. At 30 GHz, we detected 
57 sources above a limiting flux density of ~ 5 mjy. This 
is the deepest follow-up of any complete sample of sources 
detected at 15 GHz by the RT. 

At a flux density of 10 mjy, which is our estimated 



observed a much larger sample of sources down to mjy lev- 
els but selected at a much lower frequency (1.4 GHz). The 
potential danger of using low frequency selected samples is 
that there may exist a significant population of sources with 
steeply rising spectra towards high frequencies that are not 
present in the low frequency surveys. As the two surface 
density estimates are consistent, this indicates that such a 
population is not obviously present at the 10 mjy level. 

We have compared our flux density measurements with 
those from the VSA source subtractor and VLA measure- 
ments. These comparisons give confidence in our flux scale 
but reveal that a significant fraction of sources are variable 
on a timescale of a few years, some at the level of a factor 
of 2. This shows the importance of taking contemporane- 
ous measurements of discrete sources in conjunction with 
measurements of the CMB. 

We have also investigated the dependence of the spec- 
tral index distribution on flux density by comparing our 
measured spectral index distribution with that for much 
stronger sources (above 1 Jy) selected from the WMAP 22 
GHz catalogue. We conclude that the proportion of steep 
spectrum sources increases with decreasing flux density. This 
is qualitatively consistent with models of source populations, 



for example de Zotti et al. ( 2005 ) 
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• J0936+3129 - Double source in FIRST map, RT source 
coincident with western component. 

• J0941+3126 - the FIRST and NVSS maps show that 
this source has two extended components. This source is not 
used in the spectral analysis. 

• J0941+3226 - NVSS and FIRST show 15 GHz compo- 
nent is at the southern end of an extended source. 

• J1235+5317 - 15 GHz position is ~ 55 arcsec south-east 
of NVSS peak. 

• J1235+5311 - no 1.4 GHz sources are detected in the 
NVSS or FIRST maps at the RT position but weak extended 
components straddle this position in the NVSS map. 

• J1238+5249a/b - are separated by approximately 
2.5 arcmin east- west and thus for half the on-off observa- 
tions of each source the other lies within the reference beam. 
The correction, however, is very small. 

• J1535+4142 and J1535+4143a/b - lie within 3 arcmin 
of each other. Examination of the FIRST map of the region 
suggests that they may be physically related. Our 30- GHz 
flux density measurements are confused. Not used in the 
spectral analysis. 

• J1541+4114 - NVSS and FIRST maps show that 
15 GHz component is probably the central component of 
an asymmetric triple source. 



This paper has been typeset from a T[hX/ L^TjtX file prepared 
by the author. 



APPENDIX: NOTES ON INDIVIDUAL 
SOURCES 

• J0014+2852a/b - probably two components of a single 
source separated by 75 arcsec and our listed 30 GHz flux 
densities are blends of both components. Not used in the 
spectral analysis. 

• J0301+2541/42 - are separated by ^1 arcmin and are 
probably part of the same source. Their flux densities are 
blended in our 30 GHz measurements. Not used in the spec- 
tral analysis. 

• J0304+2659 - 15 GHz position is ~ 15 arcsec south of 
NVSS peak. 



